Optical transition of the 229 Th nucleus in a solid-state environment 
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We describe a novel approach to directly measure the energy of the narrow, low-lying isomeric 
state in 229 Th. Since nuclear transitions are far less sensitive to environmental conditions than 
atomic transitions, we argue that the 229 Th optical nuclear transition may be driven inside a host 
crystal with a high transition Q. This technique might also allow for the construction of a solid- 
state optical frequency reference that surpasses the precision of current optical clocks, as well as 
improved limits on the variability of fundamental constants. Based on analysis of the crystal lattice 
environment, we argue that a precision of 3 x 10~ 17 < A/// < 1 x 10 -15 after 1 s of photon 
collection may be achieved with a systematic-limited accuracy of A/// ~2x 10 _ie . Improvement 
by 10 2 — 10 3 of the constraints on the variability of several important fundamental constants also 
appears possible. 



The enormous effort towards building next-generation 
atomic clocks based on optical transitions [l|-|4[ testifies 
to the impact a more accurate time reference will have. 
Already, atomic clock experiments have provided some 
of the most stringent tests of General Relativity [j| and 
produced the tightest constraints on present day varia- 
tion of many of the fundamental constants 0] . Despite 
their successes, traditional clock experiments are cumber- 
some: the high-quality (Q) oscillator is an atomic tran- 
sition, which is sensitive to its environment. To mitigate 
environmental influences, modern clock experimenters 
routinely employ complicated interrogation schemes, e.g. 
atomic fountains or 3-dimensional optical lattice confine- 
ment. An interesting shift in paradigm is to consider 
an optical clock based on a nuclear transition. Just as 
in atomic clocks, the high-Q oscillator, i.e. the nuclear 
transition, can be addressed by laser spectroscopic tech- 
niques - as long as the transition energy is accessible with 
current laser technology - with the added benefit, due to 
screening by the atomic electrons, of reduced sensitiv- 
ity to the environment. The low-lying isomeric state in 
229 Th provides such a transition. 

The existence of a low-lying state in 229 Th, with an 
estimated natural linewidth, r„, in the range 0.1 < 
r„/27r < 10 mHz 0,0, was indirectly established over 
3 decades ago through measurements of 7-ray energies 
resulting from a-decay of 233 TJ Subsequent indirect 
measurements determined the isomeric state (I = 3/2 + ) 
to be 3.5 eV above the ground state (I = 5/2 + ) [To| ; 
however, direct searches for UV emission from this Ml 



transition were unsuccessful [111 . 1 1 21 ] . Recently, a new in- 
direct measurement has established that the transition 
energy is instead 7.6 ± 0.5 eV [Hj]. Although well into 
the vacuum ultraviolet (VUV), this energy is amenable 
to study by laser spectroscopy. 

Here, we propose a novel approach to directly measure 
the transition energy to higher precision by interrogat- 
ing 229 Th that has been doped into a VUV transpar- 



ent crystal 14, This approach has several advan- 

tages. First, in contrast to current proposals [ljl 17 1 
using laser-cooled ions that can address densities of at 
most 10 8 229 Th nuclei/cm 3 , this method results in sam- 
ples containing densities as high as 10 19 229 Th nuclei/cm 3 
(<100 /iCi/cm 3 ). The ability to address such large num- 
bers of absorbers leads to fluorescence rates sufficient for 
the use of a synchrotron light source to directly mea- 
sure the transition energy. Given that the transition 
energy is only known with ~25 nm accuracy and is in 
the VUV, this is a dramatic experimental simplification. 
Second, as is known from Mossbauer spectroscopy, nu- 
clear transitions are far less sensitive to their surround- 
ings than atomic transitions. Thus a sample of 229 Th 
atoms doped into a suitable host crystal retains a nar- 
row transition linewidth, potentially allowing for the con- 
struction of a solid-state optical frequency reference with 
a higher precision than current optical clocks. Third, 
229 Th atoms in a crystal lattice are confined to the Lamb- 
Dicke regime [IH - that is, the recoil energy is much less 
than the energy required to create a phonon in the lat- 
tice, and therefore the internal and external degrees of 
freedom are decoupled. As a result, the transition en- 
ergy will not be sensitive to recoil or first order Doppler 
effects. Nonetheless, there are still a variety of small ef- 
fects that must be considered. In what follows, we will 
show that as long as each nucleus experiences the same 
crystal lattice environment, many of these effects lead to 
constant, identical shifts rather than broadenings of the 
transition frequency. As a result, they should not affect 
the performance of a clock based on this transition. 

In the absence of external fields, the Hamiltonian for 
a 229 Th nucleus in a crystal can be written as H = 
Ho + Hps + HhfS; where Ho is the single-particle nu- 
clear Hamiltonian, e.g. the Nilsson Hamiltonian fljj ]. 
while Hps and Hhfs ar e the Hamiltonian describing 
the fine and hyperfine interactions, respectively, aris- 
ing from the crystal environment. As will be discussed 
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later, we are concerned only with ionically bonded crys- 
tals with no unpaired electron spins, therefore the fine- 
structure interaction vanishes. Thus, the hyperfine in- 
teractions are the dominant source of shifts and broad- 
enings of the nuclear transition. Using a multi-pole ex- 
pansion, the hyperfine Hamiltonian can be written as 
Hhfs = Heo + Hmi + He2 + • • • ■ These are the same 
effects that arise in Mossbauer and NMR spectroscopy, 
and as a result, each can be reliably estimated using pre- 
vious measurements. In what follows, we explain the 
estimated effect of the first few terms in the series for a 
generic lattice environment. 

The electric monopole term, Heo, describes the contact 
interaction between the nucleus and its electron cloud. 
Thus, changes in the size of either the electron cloud or 
the nucleus lead to changes in the transition energy. Since 
the electron density at a given lattice site is temperature 
dependent, the transition energy will shift with temper- 
ature, and any temperature gradients across a sample 
containing 229 Th can lead to inhomogeneous line broad- 
ening. Mossbauer spectroscopy studies indicate that this 
shift will be <100 MHz compared to free ions, and that 
the temperature dependence of the shift will be on the 
order of 10 kHz/K W 



Effect 


T-dependent Shift 


Broadening 


Monopole 


5 x 10" 12 K _1 


5 x 10" 12 K _1 


Magnetic Dipole 





<5 x KT 12 


Electric Quadrupole 


<1 x KT^K" 1 


<1 x Kr la K -1 


Second Order Doppler 


<1 x lcn^KT 1 


<1 x 10 _15 K _1 
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If the atoms reliably sub- 
stitute for a single species in a crystal the shift will be 
the same for all 229 Th nuclei. Additionally, as labora- 
tory techniques can lead to temperature stability better 
than 0.1 mK across the crystal [23|, this term contributes 
~1 Hz to the transition lincwidth. 

The magnetic dipole term, Hmi, will give rise to a shift 
and splitting of the ground and excited state energies if 
there is a non-zero magnetic field at the nucleus. Crys- 
tals that are transparent in the VUV, and thus suitable 
hosts for these studies, necessarily have large conduction 
band gaps and therefore will be free of spurious currents. 
Furthermore, in ionically bonded crystals of interest, tho- 
rium substitutes as a Th 4+ ion (isoelectronic with radon), 
so all electron spins will be paired. Therefore, the only 
magnetic fields in the crystal will be those due to the 
magnetic moments of other nuclei. The magnetic field 
produced by one nuclear magneton at the shortest typ- 
ical distance (~2.5 A) in a crystal is ^400 mG. As the 
magnetic dipole moment of the 5/2 + and 3/2 + states are 
estimated to be 0.45 ^ln and -0.08 hn respectively [3], 
this would lead to a broadening of the transition energy 
by ~400 Hz. Given the contributions of all neighbor- 
ing nuclei, this could lead to a broadening of 1-10 kHz 
and thus could be the dominant source of broadening. 
However, it appears this broadening can be reduced by 
a factor of 20 with the Mossbauer analog [HI of "magic 
angle" spinning used in NMR. 

The electric quadrupole term, He2, can also give rise 
to a shift and splitting in the transition energy if an elec- 
tric field gradient is present since both states have spin 
I > 1/2 and thus support an electric quadrupole moment. 
The electric field gradient at the nucleus depends sensi- 



TABLE I: Summary of relevant fractional shifts and broaden- 
ings for LiCAF. For broadenings, K _1 refers to the magnitude 
of the temperature variation across the sample. 



tively on the structure of the host crystal and the effect 
of this term must be considered on a case by case basis. 
Further complicating a reliable estimate of the effect is 
that only the ground state quadrupole moment of 229 Th 
is known (4.3 barns (2(|). Nevertheless, an estimate for 
the size of the effect can be obtained by considering only 
the ground state shift. Field gradients in crystals are typ- 
ically in the range of 10 16 -10 21 V/m 2 [HQ, leading to 
transition energy shifts of 1 kHz-100 MHz. However, as 
long as the 229 Th nuclei reliably substitute into a single 
lattice site, this interaction results in a constant shift of 
the transition frequency - with a temperature (T) depen- 
dence of A/// = a AT where a is the crystal coefficient 
of linear expansion. We note that He2 is minimized for 
host crystals with cubic or icosahedral symmetry [30| | . 

Because nuclei can vibrate in their lattice sites these 
interactions could lead to a broadening of the transition 
frequency if the zero-phonon transitions are mode de- 
pendent. Such dependence is likely only for the optical 
phonon modes where the 229 Th nuclei are vibrating rel- 
ative to their nearest neighbors, however, even at room 
temperature these modes are expected be frozen out [2^ | . 
Further, such an effect would lead to resolvable sidebands 
if the shift is more than the linewidth. 

Finally, the vibration of the nuclei also leads to second 
order Doppler effects which can cause a shift in the tran- 
sition energy. Using a simple harmonic oscillator model 
this effect can be estimated as A/// = 3kT/(2Mc 2 ) 
yielding a temperature dependent shift on the order of 
1 Hz/K; A/ < 0.1 mHz across the sample for the as- 
sumed temperature stability and uniformity. 

While the above analysis shows that the crystal lattice 
environment causes large, perhaps incalculable, shifts 
in the transition frequency, this modified transition fre- 
quency (in the absence of defects) is no less fundamental 
than the atomic frequencies currently used for time keep- 
ing purposes. Thus, the only effects limiting the useful- 
ness of this transition for a primary frequency standard 
are the intrinsic broadenings and environmental effects 
that shift the nuclear energies, summarized in Table U for 
LiCAF (see below). The former will impact the preci- 
sion of the frequency measurement, while the latter im- 
pacts its accuracy. We show below that even in the worst 
case scenario (10 kHz broadening due to the magnetic 
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dipolc term), the large number of absorbers allows the 
frequency to be measured with a higher precision than 
current clocks. As standard laboratory techniques lead 
to suppression of electric and magnetic fields to below 
1 mV/cm and 1 /iG respectively, the most important 
environmental effect appears to be the crystal tempera- 
ture; a 1 mK temperature difference results in a fractional 
shift of the frequency by 6xl0~ 15 . With state of the art 
thermometer accuracies of ~0.04 mK (3l|, a clock based 
on this architecture may not possess the necessary accu- 
racy to be a primary standard; however, it should display 
high short-term stability which would make it useful as 
a flywheel oscillator that could replace e.g. the hydrogen 
mascr standard. 

In addition to its potential use as a frequency refer- 
ence, the transition also displays enhanced sensitivity 
to fundamental constant variation com par ed to measure- 
ments based on electronic transitions 32], |33j . The de- 
pendence of the transition frequency on the fine-structure 
constant (a) and the ratio of the light and strange 
quark mass to the QCD length scale (X q = m q / Aqcd 
and X s = m s / Aqcd, respectively), is derived from Ho 
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as 



33J, |3J: A/ 



5 ) where c 



is a model dependent constant. Using the constants 
derived in Ref. (33j . a variation in a, X qi or X s of 1 
part in 10 1 will lead to a change of the transition fre- 
quency by ~150 Hz, ~3 kHz, and ~37 kHz, respectively. 
Hence, measurements of the transition frequency at the 
1 Hz level (see below) separated by one year could de- 
termine independent fractional variation of these funda- 
mental constants at the <10 -19 yr -1 , 10~ 20 yr" 1 , and 
10" 21 yr" 1 levels respectively; a 2-3 order of magnitude 
improvement over current limits @, H^, [3(| . 

The first steps towards these exciting applications are 
identifying suitable host crystals and determining the 
transition frequency to a greater precision. The key re- 
quirements for the host crystal are that it be reasonably 
transparent in the VUV, have a pure crystalline struc- 
ture, and chemically accept Th 4+ ions. Numerous candi- 
date hosts can be identified in the extensive family of 
fluoride crystals. Potential crystals include Na2ThFe, 
LiCAF (LiCaAlF 6 ), LiSAF (LiSrAlF 6 ), YLF (LiYF 4 ), 
and CaF2, which can be grown by various techniques, 
e.g. Czochralski, Bridgman, etc. 37, 38|. While these 



crystals have band gaps of >10eV, impurities and color 
centers can cause residual absorption and may inter- 
fere with the 229 Th optical spectroscopy. Various metal 
ions as well as oxygen-based impurities at the part-per- 
million level can lead to absorption coefficients of up to 
0.1cm -1 [3^|. To ensure crystal transparency, standard 
purification methods can be employed to reduce these im- 
purities to part-per-billion levels [401 ] . We are currently 
applying these methods to the preparation of high purity 
metal fluoride precursors for crystal growth. 

Thus far our experimental effort has focused on LiCAF 



which is transparent down to 110 nm and has all electron 
spins paired. 229 Th 4 + ions are expected to substitute at 
the Ca 2+ location in LiCAF along with two F~ intersti- 
tials for charge compensation [4l[ . Due to the octahedral 
symmetry of the Ca site, electric field gradient effects will 
be minimized. We have fabricated a Th:LiCAF test crys- 
tal using the less expensive 232 Th isotope. A density in 
excess of 10 18 232 Th/cm 3 was measured using Rutherford 
back-scattering. We have also made measurements of the 
fluorescence backgrounds which might arise from the a- 
decay of 229 Th (4.8 MeV, r = 7880 yr) by bombarding 
this crystal with a 30 nA, 4.8 MeV beam of a-particles 
at the Ion Beam Materials Laboratory of Los Alamos Na- 
tional Laboratory. A total fluorescence rate correspond- 
ing to 0.3 photons per a-particle was observed over the 
unfiltcred PMT bandwidth (Hamamatsu R7639; 115 nm 
- 230 nm). Thus, a crystal containing 10 18 229 Th nulcei 
yields a background fluorescence rate of 2 MHz. Even 
in the unrealistic case that all of this fluorescence occurs 
exactly at the transition frequency we show below that 
it is not problematic. 

To make a higher precision direct measurement of the 
transition frequency we have planned a multi-phase ex- 
periment. The first phase will be a low resolution study 
using the Advanced Light Source (ALS) at Lawrence 
Berkeley Labs. The ALS provides tunable VUV light 
(5 eV-30 eV) with a lincwidth (A) of 0.175 eV, a photon 
flux of ~10 20 photons/cm 2 -s, and a 170 /im x 50 [im 
spot size. Assuming both the light (A) and broad- 
ened transition (r) lincwidths to be Lorcntzian, then 

'fs is given 



far from saturation, the excitation rate 



j— $ p Ng where w (A) is 



b v ps ( A) -EI™ 

"* dt \2tvJ r+A 1+4( ^o^)2 

the transition frequency (wavelength), W£ is the center 
frequency of the ALS beam, and N g is the total number 
of ground state 229 Th nuclei addressed by the ALS beam. 
To locate the transition the ALS center frequency will be 
scanned across the interesting range of (7.6 ± 0.5 eV) in 
steps of 0.05 eV with 200 s of illumination followed by 
100 s of light collection at each point. Even in the lowest 
fluorescence case (T n = 2tt x 100 /iHz), we expect a min- 
imum total fluorescence rate of ~15 kHz when the ALS 
is on resonance. Assuming the a-induced background is 
the dominant noise source this rate allows for a signal 
to noise ratio (SNR) > 10:1. In this manner, the entire 
spectral region can be scanned in a few hours. Because 
the transition is both narrow and long-lived, backgrounds 
from scattered light and backgrounds from induced spu- 
rious fluorescence, which is in general relatively short- 
lived and broadband, are not expected to be problem- 
atic. Once the resonance is found, longer illumination 
will lead to sufficient fluorescence rates for a VUV spec- 
trometer measurement. We expect that the transition 
frequency can be determined to ~50 cm" 1 (0.1 nm) in 
these experiments. 

Radiation trapping will limit the total density of 



4 



229 Th to n max = 1/erw, where a is the absorp- 
tion cross section and w is the excitation beam 
waist (n max ~ 10 13 229 Th/cm 3 for the ALS measure- 
ment). However, radiation trapping can be elimi- 
nated by applying a sufficient magnetic field gradient 
(dB/dx ~ 100 G/cm) across the crystal. The concomi- 
tant line broadening is acceptable in the initial phases of 
coarsely determining the transition energy. 

After the ALS frequency measurement, we will more 
precisely measure the transition energy via laser-induced 
fluorescence. A VUV light generation system can be con- 
structed by coupling the output of a pulsed dye laser into 
a H2-Raman cell. With this system the ALS determined 
frequency range can be covered with a laser pulse energy 
of 10 /iJ and linewidth of 0.1 enr 1 [I^], allowing the 
entire range to be scanned in 6 h with SNR > 10:1. 

Following these measurements, evaluation of the tran- 
sition as an optical frequency reference can begin. The 
optimum method of precisely determining the transi- 
tion frequenc y is to probe the nuclei with a VUV fre- 
quency comb [431 ] . Locking the frequency comb to the nu- 
clear transition allows for direct comparison to standard 
clocks. Assuming the transition linewidth is 0.5-10 kHz, 
a 100 pW on-resonance comb tooth leads to a transition 
frequency resolution of 3 x 10~ 17 < A/// < 1 x 10~ 15 
with 1 s of photon collection. Longer integration times 
improve these limits, but only to the level allowed by 
temperature stability. With the best temperature stabil- 
ity demonstrated to date [23[ the resolution is limited to 
2 x 10~ 18 . This level is reached in 5 min - 5 d for the 
linewidth range considered. As mentioned above, the 
long-term accuracy is 2 x 10 -16 as set by current ther- 
mometer accuracies. 

In summary, we have proposed a novel method for the 
direct measurement of the energy of the low-lying iso- 
meric state in 229 Th. The analysis presented indicates 
that it should be possible to retain a narrow transition 
linewidth even when the 229 Th nucleus is doped into a 
crystal. The technique is limited by temperature accu- 
racy and stability which allows for a fractional frequency 
accuracy of A/// ~ 2x 10~ 16 in the long-term and a res- 
olution of Af/f ~ 2 x 10~ 18 in the short-term. The tech- 
nique looks promising for a solid-state optical frequency 
reference, as well as tests of fundamental constant vari- 
ation. The sensitive temperature dependence might also 
make it a useful thermometer. 

The ALS is supported by the U.S. DOE under Contract 
No. DE - AC0205CH11231. This work is supported by 
UCLRP award No. 09-LR-04-120497-HUDE. 
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